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Abstract. We discuss the deviations from the ideal relativistic thermal bath. These deviations are dominated from quantum 
chromodynamic (QCD) corrections in the most part of the parameter space of the Standard Model (SM) and the Minimal 
Supersymmetric Standard Model (MSSM). These effects are relevant for astrophysical precision measurements and dynamics 
of scalar tensor theories (SST). 
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THE THERMAL UNIVERSE 

The thermodynamics of any substance can be character- 
ized by its free energy density (^), or equivalently gf, 
the effective free energy number of relativistic degrees of 
freedom: 
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that is defined by normahzing to the free energy density 
of one non-interacting massless bosonic degree of free- 
dom. For vanishing chemical potential, all the main prop- 
erties of the fluid can be deduced directly from Eq. ([1}, 
such as pressure (P), energy density (p ) or their ratio { (£)): 



P 
P 

CO 







dp 






P=—T 




90 




' d\ngf 


P \ 


, + d\nT 



dlnT 



(2) 



The partition function and gf can be computed diagram- 
matically The one loop vacuum diagram provides the 
free theory result. In this case, the coefficient gf is sim- 
ply the number of bosonic (Nt,) and fermionic {Nf) rel- 
ativistic degrees of freedom: g/.free = Nb + ^ Nf, where 
the 7/8 coefficient multiplying the fermionic contribu- 
tion is due to the different value of the partition func- 
tion obtained from the Fermi-Dirac thermal distribution 
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rather than the Bose-Einstein distribution. Higher loops 
account for the interactions. In [1], it has been shown 
that the strong interaction is the dominant even at very 
high temperatures except for the case of a heavy Higgs 
boson. Unfortunately, the situation is not so simple since 
the finite temperature perturbation approach to QCD has 
a very slow convergence. In any case, the first order con- 
tribution has the right magnitude and the leading order 
correction coming from the strong interaction estimates 
the interaction effects in the thermal bath. In particular, 
for the standard model Uj]: 
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with as (m,) ~ 0.0172; whereas for the MSSM fl]: 
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Assuming for definiteness that the supersymmetric part- 
ners have masses msusY = 500 GeV, and using the run- 



ning in Eq. ^ from T — niz to T — wsusy> we have 
a3('«susY)^ 0.0153. 

We recall that the expressions (|3]l and (|4]l hold for 
temperatures greater than the masses of the particles, i.e. 
for T > 200 GeV in the SM, and for T > 500 GeV in 
the MSSM with our assumption. At low temperatures, 
the effect is stronger in the SM than in the MSSM case. 
This is because 0:3 runs faster in the SM. However, since 
a3 decreases less in the MSSM, the departure from the 
noninteracting value becomes greater for the MSSM as 
the temperature increases. 



SCALAR-TENSOR MODELS 

The deviation from the ideal gas is particularly impor- 
tant in scalar-tensor theories (STT) of gravity, since it is 
not only a correction but also the main contribution to 
the scalar evolution. Many different types of new scalar 
fields are predicted different extensions of the standard 
model of particles and general relativity [3,4], but a stan- 
dard STT can be defined by the introduction of a spin-0 
scalar mediator in addition to the standard spin-2 gravi- 
ton. Well known examples of such fields are the Jordan- 
Fierz-Brans-Dicke (JFBD) scalar which is introduced as 
an extension to general relativity (GR), the graviscalar 
(or radion) in extra dimensional models, or modili fields 
(as the dilaton) in the context of string theory. 

The general action of these models can be written as: 
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where the first term is the action for the spin-2 graviton 
and the scalar field, while the second term is the action 
for matter This expression defines the theory in the Jor- 
dan frame, in which the standard expression of the met- 
ric gi^v is used in the action for the matter fields (where 
by "matter" we generally denote any field in the theory 
apart from and gjiv)- However, the gravitational inter- 
action between matter fields is modified, since the func- 
tion A{(p) multiplies the Ricci scalar Such a theory de- 
pends on two arbitrary functions. For the sake of sim- 
plicity, we will assume that G*, the bare gravitational 
constant, is approximately the gravitational constant that 
would be measured in a Cavendish-type of experiment. 
Under a conformal transformation 
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and for a vanishing potential (U = 0), the action of the 
system becomes 



S = 
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where the scalar field 0* is defined by 
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and A (0*) is short for A (0 (0*)) . 

The expression ^ is the action of the system in the 
Einstein frame. It is characterized by a standard ac- 
tion for the spin-2 graviton; however, the combination 
A^ {(j>t)g*^v, rather than the metric itself, is used in the 
action for matter While the two expressions are equiva- 
lent, the Einstein frame is more often used to study the 
cosmology of the system (since the resulting cosmolog- 
ical equations are the standard ones), while the Jordan 
frame is more often used to study particle physics pro- 
cesses (since the physical lengths and masses are con- 
stant in this frame). Moreover, the Einstein frame is built 
in such a way that the kinetic term of the spin-2 and spin- 
mediators is diagonal so that the Cauchy problem is 
well-paused in this frame. 

The cosmological evolution of the scalar mode is de- 
termined by the standard equation of motion for a scalar 
field with an extra source term proportional to the trace 
of the energy momentum tensor Assuming a perfect re- 
latisvistic fluid form for the energy momentum tensor, its 
trace is simply, p — 3p, and vanishes in a radiation dom- 
inated epoch when the equation of state is characterized 
by p = p/3. While this is certainly a good approximation 
deep in the radiation era. A non-vanishing source term 
also arises when one includes the contribution from the 
trace anomaly. Although this is generally small, typically 
w~l/3 = & (lO""* - 10"^) li], there can be situations 
in which this contribution has a significant impact on the 
cosmology of these models. In fact, not only can it af- 
fect the evolution of the scalar field, but in some cases, it 
can lead to a brief early phase of contraction followed by 
a non- singular bounce in the cosmological scale factor 
In these cases, one can also derive a maximum tempera- 
ture of the universe which may provide a solution to the 
gravitino problem or other unwanted relics [IJ. 



FINAL REMARKS 

We have analyzed the consequences of deviations from 
the ideal relativistic thermal bath for the dynamics of 
scalar-tensor models. However, interaction effects are 
also important in the precision era of astrophysical ob- 
servations that we have already entered. This fact can be 
illustrated in the case of a thermal relic. Although there 
are other possibilities [5], Dark Matter (DM) is usually 
assumed to be in the form of stable Weakly-interacting 
massive particles (WIMPs) that naturally freeze-out with 



the right thermal abundance. WIMPs emerge in differ- 
ent well-motivated particle physics scenarios as in R- 
parity conserving supersymmetry (SUSY) models llil 
[tI], universal extra dimensions (UED) [8, 9], or brane- 

worids iii3,[ni[il. 

The present uncertainty on the total DM density is 3% 
ifTsll . and it can be improved below 1 % by analyzing data 
from the Planck Surveyor. However, in order to associate 
the observed precision with the thermal relic density of a 
particular WIMP model, it is necessary to understand the 
thermal bath up to the same level. Unfortunately, the un- 
certainties are bigger than 1 % percent even at high tem- 
peratures. Therefore, non-perturbative analises, as lattice 
studies, are necessary. Alternatively, if we are able to un- 
derstand the nature of DM from other experiments as the 
new generation of colliders lfl4l . we may have the oppor- 
tunity to improve our knowledge about QCD by making 
precise astrophysical observations. 
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